Different length-scales for order parameters in two-gap superconductors: 
the extended Ginzburg-Landau theory 
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Using the Ginzburg-Landau theory extended to the next-to-leading order we determine numeri- 
cally the healing lengths of the two order parameters at the two-gap superconductor/normal metal 
interface. We demonstrate on several examples that those can be significantly different even in the 
strict domain of applicability of the Ginzburg-Landau theory. This justifies the use of this theory 
to describe relevant physics of two-gap superconductors, distinguishing them from their single-gap 
counterparts. The calculational degree of complexity increases only slightly with respect to the 
conventional Ginzburg-Landau expansion, thus the extended Ginzburg-Landau model remains nu- 
merically far less demanding compared to the full microscopic approaches. 
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Over the past half-century, the Ginzburg-Landau (GL) 
theory^ has proven to be a very helpful tool in studies 
of superconductors, but also other systems in and out 
of the condensed matter physics. By its construction, 
the formalism is only justified near the critical temper- 
ature T c , but it typically produces qualitatively correct 
results even far below the T c . Recently an exception was 
found to this unwritten rule, when it was shown that 
the standard formulation of the GL theory is insufficient 
for adequate description of two-band (or multi-band) su- 
perconductors because it predicts the same spatial vari- 
ation of the condensates in all bands 2 (within its range 
of applicability i.e. close to T c ). This renders it unable 
to connect to the theoretical results obtained by using 
Usadel 3 -^ or Eilenberger— equations in the broader tem- 
perature range, which unambiguously show presence of 
two separate length scales for two gaps. Experimentally, 
the evidence for different coherence lengths in a and 7r 
bands of MgB2 was previously found by the direct vortex 
imagingS. and in pronounced features in the behaviour of 
the thermal conductivity^, specific heat£ and flux-flow 
resistivity^ as a function of applied magnetic field. 

In order to capture the important physics of different 
length scales, one needs to extend the GL formalism as 
realized in Ref. [ToL where the two order parameters are 
calculated up to order r 3 / 2 in the small deviation from 
the critical temperature r = 1 — T/T c , instead of the 
standard r 1 / 2 as used in Ref. 0- In the latter article, 
the authors mention that the extra terms in the next-to- 
leading order are by construction small corrections, not 
significant enough to alter the single coherence length 
controlling the spatial distribution of both condensates. 
First, we argue that the effects of higher-order correc- 
tions can be significant since the above argument about 
small corrections applies only to the order parameter - 
not necessarily to its spatial profile i.e. healing lengths 
and other physical quantities. Second, we argue that any 
(even small) difference between the characteristic length 
scales of two condensates is of fundamental importance, 
since it may lead to other novel phenomena related to 
the competition of length scales. Note that here we do 



not enter the recent debate about sufficient discrepancy 
of length scales to provide 'type-1.5' superconductivity^. 
Instead, we complement that discussion, by exactly quan- 
tifying the difference in length scales in the domain of the 
extended GL theory. 

The fact that the difference of characteristic length 
scales of two Cooper pair condensates exists and can 
be significant even in the strict Ginzburg-Landau do- 
main is of great practical importance, since the calcu- 
lations based on the Ginzburg-Landau theory are typi- 
cally far less computationally demanding than the calcu- 
lations based on full microscopic theories (Bogoliubov-de 
Gennes, Gor'kov, Usadel or Eilenberger equations). In 
microscopic formalisms, one usually has to make clever 
approximations and make compromises in the calcula- 
tion procedure. As a consequence, even though these ap- 
proaches are valid in the whole temperature range from 
absolute zero to T c , they are limited to the simple sys- 
tems such as a single vortex or other highly symmet- 
ric or effectively one-dimensional cases. On the con- 
trary, the Ginzburg-Landau equations have a much sim- 
pler structure and therefore allow for studying of highly 
non-trivial situations such as vortex lattice statics and 
dynamics, current driven systems, interaction with pin- 
ning and fluxonic devices. Since the above is well estab- 
lished in the standard GL formalism, we emphasize here 
that the calculations become only a fraction more com- 
plicated in the extended Ginzburg-Landau formalism - 
while it does contain more equations (for the order pa- 
rameters to the leading, and next-to- leading order), the 
coupling of those equations is realized through the coef- 
ficients and not the calculated variables - contrary to the 
standard set of two-band GL equations^. 

Theoretical formalism Following Ref. [ToL we employ 
the extended GL formalism for the order parameters Aj 

evaluated up to order t 3 / 2 by taking A 3 -(x) = A^(x) + 



A^(x), where A) UJ (x) oc t 1 / 2 and A^(x) oc r 3 / 2 , with 
j = 1,2 indexing two coupled condensates in a two-band 
superconductor. In the absence of applied magnetic field, 
the order parameters can be taken real. The extended GL 



2 



equations then read (at zero magnetic field) 
aAf +/3 i [Af] 3 -i07 2 Af = 0, 



analysis shows that tAj ~ Ay. Then, by construc- 
tion, Aj. > A*- and the trivial inequality r < 1 holds. 
However, this is a necessary but not sufficient condition 
justifying Eqs. @, as one also needs to keep in mind the 
justification for the use of the gradient expansion in the 
derivation of Eqs. ©• This expansion requires the GL 
coherence length £ to be larger than the band-dependent 
correlation length Q — hvj/(2irT) that controls the spa- 
tial variations of the relevant kernels in the integral ex- 
pansion of the anomalous (Gor'kov) Green function in 
powers of the band-dependent order parameters (in the 
clean limit), i.e., 
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where a = — N(0)t (ni_4 2 + ^2-4i)/Ai 2 , K = 
N{0)h 2 Wl{ ni vlA 2 + n 2 t- 2 2 A)/(6A 12 ), /3 1(2) 
7V(0)T^ 3 2 (n 1(2) ^ 2(1) /A 12 + n 2 (i)-4? (2) /A? 2 ), with 

Ai = A 22 — ijniA and A 2 = An — nn 2 A. Here 
A = \n(2e r hw D /irT c ), with Euler constant T= 0.577, 
and rj denotes the determinant of the interaction matrix 
Ay = N(0)gij, with g^j the coupling constant, jV(O) the 
total density of states (DOS) and rijN(0) the band- 
dependent DOS. In addition, the coefficients feature 



Fermi velocities Vj of both bands, W- 
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jk, with 



£(. . .) the Riemann zeta- function, and terms F and 
Fj which are given in complete form in Ref. [Tol Since 
Eq. (QJi) is completely equivalent to the single-gap GL 
equation and the bands share the same critical tempera- 
ture T c , we have either both bands normal or both bands 
superconducting. In the first case, Eq. (QJ>) allows 

only A*; 1 ' = as a solution and thus superconductivity 
cannot be restored by corrections in the extended GL 
model. In the following we will consider that both 
bands are superconducting and their bulk amplitudes 

Wj- 

Af } = WjAf 



y/—a/fij are real. We rescale the equations using 

and x = £x, where £ = \J —Kj a is the 
length unit common for both condensates, and obtain 
(tildes omitted) 



A 
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[a(°)]3 + v 2 a( 0) = o, 



A( 1) (f-3[Af ) ] 2 )+V 2 A 



, F{Af)+F 2 {Af) 
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The first equation is the same for both bands, while 
the second one differs through the terms on the right 
side. This is exactly the cause for the emergence of two 
different characteristic lengths in the two bands. From 
Eq. ([2]) it is also directly apparent why this does not 
happen in the order parameters of order t 1 / 2 , i.e. in 
Eq. (|2^), but only in the higher order considerations as 
in Eq. (0d). To evaluate the coefficients entering Eq. @ 
one needs in principle to specify the coupling constants 
An , A 22 and Ai 2 , the partial density of states in one band, 
e.g. ni (since n\ + n 2 = I this determines n 2 as well), 
and the ratio of the Fermi velocities t>i/v 2 . All other 
parameters enter units of scaling, and therefore have no 
impact on any physical effects. 

The Ginzburg- Landau domain Before discussing the 
numerical results of Eqs. @, it is very important to get 
a feeling about the relevant values of r, i.e. the domain of 
applicability of the extended GL model. First of all, our 



Ci,C 2 <£. 



(3) 



It is important to note here that for a two-band supercon- 
ductor Q is not necessarily close to the band-dependent 
generalization of the Pippard length, i.e., hVj /(7tAj)|t=o> 
often used as an estimate of the coherence lengths of two 
bands 2, ii. Using definitions of a, K and f3j, we obtain 
from Eq. ([3]) the following estimate for the GL domain 
t S, t* with 
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l + n 2 Ai/{n 1 A 2 ) 



(vi < v 2 ). (4) 



If vi > V2, the replacement lf)2 should be made. The 
Eq. (Ql shows that the extent of the Ginzburg-Landau 
domain for specific multi-gap superconductor depends 
on its material parameters, as was hrst pointed out in 
Ref. S The Eq. © is most sensitive to the ratio of the 
Fermi velocities vi/v 2 . In particular, when V\jv 2 = f, 
the above inequality simply implies r < 0.38, i.e., the 
same as in the one-band case, regardless of the par- 
ticular values of the other parameters, i.e., n 2 /n\ and 
A 2 /Ai (where A 2 /Ai = [A^ 0) /A 2 0) ] 2 ). For m w n 2 and 
A\ w A 2 , we typically obtain 0.27 < t* < 0.38, thus 
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FIG. 1: The estimated temperature range r < r* of the valid- 
ity of the GL theory based on Eq. . The gradient expansion 
made in the derivation of the GL equations (both standard 
and extended) is fully justified below the plotted 3D-surface. 
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the temperature domain of GL theory is still very large. 
The GL domain shrinks significantly only when the ra- 
tio n 2 Ai/(niA 2 ) acquires either extremely large or very 
small values, which, e.g., can occur when the two en- 
ergy gaps differ more than an order of magnitude and at 
the same time v\ is very different from v 2 . For example, 
for the magnesium diboride parameters (as elaborated 
further below) one finds from Eq. © r* = 0.08 - 0.32 
depending on the relevant ratio of the Fermi velocities 
(i.e. on the direction considered, because a band of mag- 
nesium diboride is highly anisotropic), which is an order 
of magnitude larger than the estimate r* pa 0.02 given 
in Fig. 1 of Rcf. 4. The threshold temperature for the 
applicability of the GL theory (r*) for a broader range 
of parameters (i.e. other materials) is plotted in Fig. [TJ 
Different length scales in selected borides, pnictides, 
and nanothin Pb films To evaluate and compare the 
spatial distribution of the two order parameters in a two- 
band superconductor, we study the simple case of a su- 
perconductor/normal metal (S/N) interface in the ab- 
sence of any applied magnetic field. It then suffices to 
consider the one-dimensional version of Eqs. ©, along 
the x-axis perpendicular to the S/N interface. The ap- 



propriate boundary conditions are A 



(fc=0,l), 



0, 



and \7Af\x 



(x = 0) 

— > do) = 0, where thus x = lies at 
the interface, with superconductivity fully suppressed 
there. The second boundary condition ensures unper- 
turbed two-gap superconductivity away from the S/N 
interface. In Fig. [2fa) we show the results for Aj(x) 
in MgB2 (normalized to bulk values, at temperature 
T = 0.95T C , and with other parameters taken from Refs. 

Hand EH), 

obtained within standard^ 2 -, extended^, and 
reduced 2 - GL formalism. In the case of MgB2, the ex- 
tended GL model clearly gives two length scales for two 
order parameters Aj = i t 2, smaller than those obtained in 
the standard GL theory (with incomplete higher-order 
terms). 

For correct comparison one must define a measure for 
the spatial variations of Aj. From the numerical so- 
lution of Eq. (J5Ji) we found that at the characteristic 

distance £ the single-gap-like order parameter A^ in- 
creases from zero to 0.6089 of its bulk value Ajq. There- 
fore we define the healing lengths £j for the two con- 
densates in a two-gap superconductor using the criterion 
Ajfe) = 0.6089A j0 , as depicted in Fig. Hfb). The dif- 
ference between the characteristic length scales in two 
bands is clear already from their definition in Fig. [2(b), 
but we emphasize this point in Figs. [2](c,d), where the 
ratio of the healing lengths of the two order parameters 
is plotted as a function of temperature [Fig. [2(c)] and 
ratio of Fermi velocities vi/v 2 [Fig. [21(d)] . At the lowest 
temperature shown in (c), T = 0.92T C , the disparity be- 
tween the healing lengths is already over 25%, and the 
difference increases as the ratio of the Fermi velocities is 
taken smaller [see (d)]. 

In the analysis of the above phenomenon we noticed 
that only the coefficient S in front of the term V 2 A^-°' in 
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FIG. 2: (a) The spatial profile of the order parameters at the 
S/N interface at T = 0.95T C , for microscopic parameters of 
MgB 2 = vi/vi = 0.255^; An = 1.88, A 22 = 0.5, Ai 2 = 

0.21, m — 0.43^), compared in three versions of the two- 
band GL theory. The ratio of the healing lengths of the two 
coupled condensates £1/^2 calculated in extended GL theory 
[with £1,2 determined as illustrated in (b)] is shown in (c) as 
a function of temperature for fixed V1/V2 = 0.255, and in (d) 
as a function of the ratio V1/V2 at T — 0.95T C . 



F(Aj) depends both on r and Vi/v2- 



8 = 



h 2 N(0)Wj' 
6A12 



E v%2n 3 X n -T 1 n 1 n 2 {l + 2A)). (5) 
3=1,2 



Due to this special form, for particular values of nj and 
Xij the term in brackets can be positive for one band, 
but negative for the other. In such a case it is possible 
to change the sign of S term by varying the ratio of the 
Fermi velocities 1)1/1)2 around its threshold value 



2n 2 X 22 - ■qn 1 n 2 {l + 2A) 
2niAn — rjnin 2 (l + 2A) 



1/2 



(6) 



which equals 0.46 for the parameters of Fig. [5] Therefore, 
two cases with v\ jv 2 significantly larger and smaller than 
0.46 will show very different behavior with respect to the 
disparity of the healing lengths of the two condensates, 
as is visible in Fig. GJd). For v\jv 2 rs 0.46, the term with 
coefficient S has no influence and other terms determine 
the spatial behavior of the order parameters. 

Besides the difference in the healing lengths of the 
two condensates, the extended GL model provides in- 
sight also in a more accurate temperature dependence 
of the order parameters. Namely, standard GL theory 
leaves one with simply A oc r 1 / 2 , which is shown in 
Fig. [3] to be clearly inadequate for our calculations. 
Specifically, in Fig. [3] we show two healing lengths as 
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FIG. 3: The healing lengths of the two order parameters in 
MgB2 [Fig. Ufa)] scaled to the standard GL length scale £ as a 
function of temperature. The lines show the fitting functions 
of the data. 
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FIG. 4: The temperature dependence of the ratio of the heal- 
ing lengths in LiFeAs, for parameters from Ref. Il6| and 

Ref. Hi 



a function of temperature, both already normalized to 
the temperature-dependent length £ oc t~ 1//2 . It is di- 
rectly obvious that we are left with a non-constant value, 
and we fitted the general temperature-dependence of the 
residue to a quadratic form = 1 + At + Bt 2 , as 
shown in Fig. [3pi 

In what follows, we extend our analysis to materials 
other than MgB 2 . We first address the recently dis- 



covered iron-pnictides, more specifically LiFeAs. Inter- 
estingly enough, several existing experimental works on 
this material provide us with radically different micro- 
scopic parameters obtained from the fit of the superfluid 
density data in the self-consistent 7-mode l 15 i 16 and from 
two-band fitting of the upper critical fiel d 17 ' 18 . In Fig. 
IU we show the GL-calculated ratio of the healing lengths 
of two bands in LiFeAs, for the parameters taken from 
latter two references. Although different, both curves in 
Fig. [4] clearly show a discrepancy between the healing 
lengths, up to 20% at temperatures above 0.8T C . 

The present study and its conclusions are relevant not 
only for bulk multi-band materials, but also for nanoscale 
superconductors, e.g. single-crystalline nanofilms^, 
where the multi-band structure appears due to quantum 
confinement. In particular, the Pb(lll) nanofilms with 
thickness 4 and 5 monolayers are two-band superconduc- 
tors due to the presence of only two perpendicular single- 
electron levels below the Fermi energy. Our preliminary 
analysis of such systems shows that the difference of two 
healing lengths therein is below 5% in the GL domain. 
However, in quantum confined systems, either nanofilms 
or pancake fermionic condensates, the exact position of 
the bands with respect to the Fermi level (and therefore 
the vi /i>2 ratio) can be tuned by confinement^ and with 
that the discrepancy between the corresponding healing 
lengths can be significantly enlarged. 

Conclusions To summarize, we demonstrated on ex- 
amples of several multi-gap superconducting materials 
that characteristic length scales of coupled condensates in 
two- or multi-band samples can be significantly different 
from each other in the domain of the extended Ginzburg- 
Landau theory. This makes our model an excellent tool 
for further studies of two-band systems, where one ex- 
pects a plethora of novel physical effects emerging from 
the competition of different length scales. 
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dation (FWO-V1), the Belgian Science Policy (IAP), and 
the ESF-INSTANS network. 
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